Langevin description of critical phenomena with two symmetric absorbing states 
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On the basis of general considerations, we propose a Langevin equation accounting for critical 
phenomena occurring in the presence of two symmetric absorbing states. We study its phase diagram 
by mean-field arguments and direct numerical integration in physical dimensions. Our findings fully 
account for and clarify the intricate picture known so far from the aggregation of partial results 
obtained with microscopic models. We argue that the direct transition from disorder to one of two 
absorbing states is best described as a (generalized) voter critical point and show that it can be 
split into an Ising and a directed percolation transitions in dimensions larger than one. 
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The classification of equilibrium phase transitions into 
universality classes by just identifying their relevant in- 
gredients (i.e. symmetries, conservation laws, and di- 
mensionalities) constitutes one of the most remarkable 
achievements of modern statistical mechanics. Ginzburg- 
Landau- Wilson free-energy functionals, either in their 
static Q or dynamic versions (usually written in terms 
of Langevin equations, as the Model A of Hohenberg 
and Halperin describing the kinetic Ising class |2j) pro- 
vide a compact and systematic theoretical framework to 
represent universality classes: Being continuous (coarse- 
grained) theories, they are thus susceptible to analytical 
studies by using the tools of statistical field theory and 
the renormalization group. 

Out of equilibrium, the situation is far from being as 
satisfactory. In spite of evidence of universality, the rel- 
evant ingredients for classification are often not known, 
and continuous descriptions in terms of Langevin equa- 
tions or dynamical generating functionals are mostly 
lacking. For instance, within the prototypical case of ab- 
sorbing phase transitions, where the ordered "absorbing" 
states are devoid of fluctuations allowing the return to the 
disordered "active" ones, the directed percolation (DP) 
class is prominent and very robust 0, 0] . Loosely defined 
as the class of all phase transitions into a single effective 
absorbing state without extra symmetries or conserva- 
tion laws, it is represented by a Langevin equation which 
can be renormalized satisfactorily 0, [|| . But such a con- 
tinuous description or even a controlled renormalisation 
procedure is lacking for the also rather well-established 
class of phase transitions into one of two symmetric ab- 
sorbing states, despite some thoughtful attempts |(J. 

In this Letter, we propose, on the basis of general con- 
siderations, a Langevin equation accounting for critical 
phenomena with two (Z2-) symmetric absorbing states. 
Pending a renormalisation group approach, we study the 
phase diagram of this equation by mean-field arguments 
and direct numerical integration and show that it fully 
accounts for the rather intricate picture known so far 



from the aggregation of partial results obtained with mi- 
croscopic models, a situation which we briefly recall now 
before proceeding with our findings. 

The lack of consensus about the characterization of 
phase transitions into two symmetric absorbing states 
is reflected by the different names given to this class in 
the literature 0, 0, II E3 ■ Sometimes simply called 
DP2 (marking the existence of two absorbing states), or, 
more accurately, "directed Ising" (referring to both the 
Ising /^-symmetry and the presence of absorbing states), 
it is most often called "parity-conserving" because of its 
usual interface representation where, in one space dimen- 
sion, diffusing particles A undergo the reactions A — ► 3A, 
2 A — ► 0. (These branching and annihilating random 
walks with an even number of offsprings stand for inter- 
faces between domains of the + and — absorbing states, 
see Fig.^.). However, there now exists ample evidence 
that the conservation of the parity (of the number of in- 
terfaces or particles) is not the relevant ingredient 0, |j| . 
Moreover, this particle representation just gives rise to 
trivial phase transitions (at zero branching rate and with 
mean- field exponents) in higher space dimensions d > 2. 
Rather, as was briefly hinted at in 0] and suggested in 
[l2|. we endorse the viewpoint that this type of critical 
phenomenon, where interfaces between the two symmet- 
ric absorbing states branch and annihilate, is best de- 
scribed as the (generalized) voter class in the sense of 
[l^ . Recall that in the usual voter model 0] , randomly 
chosen Ising spins take the value of one of their randomly 
chosen neighbors: then only interfaces (-1 — pairs) evolve, 
with -| > H — h or with equal probability i. In di- 
mension d — 2, this model is critical and at its upper crit- 
ical dimension. It is characterized by a marginal ordering 
process during which the density of interfaces (H — pairs) 
decays like I /hit In contrast, this simple rule is not 
critical in other dimensions: for d = 1, it coincides with 
the annihilation process 2 A — >■ 0, while in d = 3 it leads 
to a disordered phase. One possible generalisation of this 
"classical" voter rule preserving the Z2 symmetry is to al- 
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FIG. 1: (Color online) (a) Schematic correspondence between 
domain and interface representation of one-dimensional gener- 
alized voter dynamics. Here an interface separating a + from 

a — domain branches (A — > 3 A, swap of two spins H > — h), 

creating two new domains; later, two of the three interfaces 
annihilate (2 A — > 0), suppressing a + domain. The two sym- 
metric absorbing states impose parity conservation of inter- 
faces, (b) variation of the deterministic potential with a for 
6=1 (from top to bottom a = 0, 0.5 and 1). (c) same as (b) 
but for b = —1 and a = —1, —1/3, 0. 

low for spin swaps H > — h [which in d = 1 amounts to 

the branching reaction A — > 3A, see Fig.^a)]. It is sim- 
ple to realize that this generalisation (and others) allows 
for tuning the model at criticality in any dimension, that 
in d = 2 the critical properties of the original model are 
preserved, and that in d = 1 the critical point is nothing 
but of the DP2/directed Ising/parity-conserving class. It 
is thus both meaningful and useful to denote this class, 
in any dimension, as the generalized voter (GV) class. 

Extrapolatingfrom their numerical results in d — 2, 
Dornic et al. 13] conjectured that transitions with Zi 
symmetry and no bulk fluctuations (i.e. with two sym- 
metric absorbing states) should all display GV critical 
points. In a recent paper, though, Droz, Ferreira, and 
Lipowski 0| somewhat challenged this picture by show- 
ing that such a transition, in some versions of two- 
dimensional generalized voter models, may not be di- 
rect, but split into a first Ising-class, symmetry-breaking 
transition followed later by a DP-class transition to the 
absorbing state chosen after the previous Ising critical 
point. One can thus legitimately wonder whether the di- 
rect GV class transitions observed by Dornic et al. — as 
well as Droz et al. — exist at all (at codimension-1 man- 
ifolds of parameter space) or whether they are just the 
artifact of close-by Ising and DP transitions, coinciding 
only at special points like the classical voter model. 

To summarize, in two space dimensions, the question 
of the possibility of the merging of an Ising and a DP 
line into a full GV line remains open, while in one space- 
dimension the relevance of parity conservation is still de- 
bated. Below, we address both of these points and clar- 
ify the nature of all phase transitions in the presence of 



two symmetric absorbing states via the introduction of a 
unique, well-behaved, Langevin equation for this general 
problem. 

Our proposal is by no means unique, but it is con- 
strained by general guidelines: The equation has to be 
symmetric under reversal of the field (<p — > — 0), which 
takes values between two absorbing barriers, set, without 
loss of generality, at ±1 ((f) € [—1,1]). Because, in two 
dimensions, the transition can be split into an Ising and a 
DP point, each of the absorbing barriers must be similar 
to those of the Langevin equation for DP, i.e. the square 
root of the distance to each barrier must appear as a mul- 
tiplicative factor of the noise. This is also corroborated 
by the fact that the Langevin equation proposed once for 
the classical (integrable) voter model 

d t <t> = DV 2 <p + a^/l-(f> 2 n (1) 

where r\ is a Gaussian noise delta-correlated in space and 
time, was recently shown to behave as expected (i.e. log- 
arithmic decay of the density of interfaces) |l7lll£j ]. In or- 
der to represent the possibility of Ising-like spontaneous 
symmetry breaking, we need to add a minimal number 
of polynomial terms with odd powers of 4>. (At least two 
free parameters are needed to describe for the splitting 
scenario uncovered by Droz et al.). We are then almost 
ineluctably led to the following equation: 

d t <j> = (a(f> - 60 3 )(1 - cf 2 ) + DV 2 (t) + a^/l-(j) 2 r) (2) 

Note that removing the 1 — 4> 2 factors, both in the de- 
terministic force and in the noise amplitude, leads to the 
Model A for the Ising class 0- Let us now describe the 
different possible regimes of Eq. in the (a, b) param- 
eter plane, a natural choice since for a = b = one 
recovers the voter equation 

We start with a discussion at the mean-field level, i.e. 
reducing Eq. to its first term, rewritten as —V'(<f>), 
with the "potential" V(4>) = -f 2 + - f^ 6 . 

b > 0: separate Ising and DP transitions. For a < 0, 
<f) = is locally stable, while it is unstable for a > 
(Fig. QJd) . The b(j) 3 term enforces stability as in Model 
A (even if the absorbing barriers are removed). At a — 
0, where the local stability around tj> = changes, the 
symmetry is broken, and we expect an Ising transition 
in the full problem |l9j. Increasing a > 0, the minima of 
the potential move progressively closer to the absorbing 
barriers and, for a — b, a collapse onto the absorbing 
barrier selected by the previous spontaneous symmetry 
breaking takes place. This second transition should be 
in the DP class once fluctuations are incorporated. 

b < 0: unique GV transition. If b = the potential is 
Gaussian around the origin, which is a stable extremum 
if a < 0, and unstable otherwise. The transition is at 
a = 0, but there is no c/> 4 term in the potential forcing it 
to be continuous: the location of the potential minimum 
changes abruptly from </> = Oto0 = ±l. This time the 
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symmetry breaking occurs simultaneously with a fall into 
one of the absorbing states. The critical point should 
be that of the voter model once fluctuations are taken 
into account. For b < 0, the </> 4 term is present in the 
potential, but it is not stabilizing, and it does not lead 
to a continuous transition. For a < the origin is locally 
stable, and there are also extrema at the barrier, and 
additional maxima at ±a/ a/b that may or may not lie 
in the interval [—1,1]. As a approaches the extrema 
move closer to the origin, and the minima at the barriers 
deepen [Fig.^c)]. At some point, the stability is globally 
changed and we expect to have a situation similar to the 
one for 6 = 0, i.e. a unique GV transition. 

The naive mean-field parameter diagram of Eq. @ 
thus consists of a symmetry-breaking line along the b- 
axis, joined by a DP line a — b at the origin. For b < 0, 
the two lines merge, and a unique GV transition occurs, 
while for b > the Ising-like symmetry breaking does 
not lead directly to one of the absorbing states, a sit- 
uation similar to that uncovered by Droz et ai. More 
elaborated, self-consistent mean-field approaches lead to 
similar results, albeit with the symmetry-breaking line 
not being along the o-axis anymore. 

In order to go beyond mean-field and elucidate the in- 
fluence of fluctuations in the phenomenology of equation 
(j2J), we integrate it numerically using the approach de- 
tailed in [la, l2i| • This method is designed to circumvent 
the numerical difficulties associated with the presence 
a singular square-root noise near an absorbing state in 
Langevin equations. It consists in separating the integra- 
tion of the deterministic terms from the stochastic piece, 
the latter being performed by sampling exactly the condi- 
tional probability distribution function (p.d.f.) solution 
of the associated (forward) Fokker-Planck equation. This 
sampled value is next used to evolve the remaining de- 
terministic part. Here, the Fokker-Planck equation for 
^ — cry/l — (j) 2 rj(t) can be solved through an eigenfunc- 
tion expansion, leading to a rather complicated p.d.f., 
with a continuous part and two delta peaks at the barri- 
ers 4> — ±1 j^jj. However, as remarked in [l^, discarding 
the (exponentially suppressed) influence of the farthest 
absorbing state, one can treat the noise term as two in- 
dependent DP barriers, and apply the existing efficient 
procedure for the DP noise. Thereby, using time-mesh 
At = 0.1, space- mesh Ax = 1, and the parameter values 
D = 0.5 and a 2 = 0.8, Eq.(J2J) can be faithfully integrated. 

We first present our results obtained in two dimen- 
sions, which agree qualitatively with the phase diagram 
predicted from mean-field arguments: 

For b larger than some b* (6* « 0.50 with our choice 
of parameters), two distinct transitions are encountered 
upon increasing a and they merge linearly as b — ► b* 
[Fig- Eta) andEtb)]. At low a, any initial condition leads 
to a disordered state ((0) = 0). For ai S i ng < a < a^p, the 
steady state, reached after some phase ordering transient, 
has a non-zero magnetization to = (</>(r,t)) r but is still 
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FIG. 2: (Color online) Results from simulations in d = 2. (a) 
Phase diagram in the (a, b) plane (error bars are smaller than 
the symbols size), (b) Steady-state magnetization (circles) 
and density of interfaces (squares) vs a for 6 = 1. (c) 1/p vs 
\nt at b = —0.2 for various values of a around ogv — —0.115 
(middle curve); the dashed line is a linar fit. (d) Time-decay 
of 1 — \m\ for a values around aop ~ 1.6551 (6 = 3); at 
criticality, 1 — \rn\ ~ t e with 9 ~ #dp ~ 0.45 (dashed line). 



fluctuating (0 < |m| < 1) and the density of interfaces 
p=l — (4>(r,t)(j)(r + e v ,t)) (where represents any of 
the unit vectors of the underlying square lattice) is finite. 
For a > odp, ordering is complete (m = ±1, p = 0). We 
have checked that the symmetry-breaking transition oc- 
curring at ai s i n g is in the Ising universality class, both by 
steady-state finite-size scaling analysis, and by measur- 
ing the decay of the time auto-correlation function from 
disordered initial conditions. For instance, the curves for 
the so-called Binder cumulant at different system sizes 
all cross each other around the universal value U* ~ 0.61 
(not shown). We have also checked that the fall into one 
of the absorbing states is a DP-class phase transition, as 
e.g. testified by the algebraic decay in time of the activity 
in a large system after a critical quench [Fig.Etd)]- 

For b < b* , a unique transition is observed at a = agv, 
across which the steady-state magnetization jumps from 
zero to ±1. On the other hand, the density of interfaces 
goes continuously to zero as a — * acv from below (not 
shown). At a = acv, the logarithmic time decay of p is 
one of the hallmarks of the GV class [Fig.Et )]- 

In one space dimension, where general arguments ex- 
clude the existence of an Ising transition, we expect a 
unique, continuous direct transition from a disordered 
phase (m = 0, p > 0) to one of the two absorbing states 
(to = ±l,p = 0). It is therefore near at hand to sur- 
mise that the ensuing critical point should be character- 
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FIG. 3: (Color online) Results from simulations in d — 1. (a) 
space-time plot of <f> at criticality a = —0.1255,6 = —0.25, 
system size 512, time running downward for 10 4 . (b) time 
decay of the interface density (p) for b = —0.25 and, from 
bottom up: a = -0.12, -0.124, -0.125, a G v ^ -0.1255, - 
0.126, -0.127; in the absorbing phase (a = -0.12), p ~ l/Vt; 
inset: data around ogv multiplied by i 286 (system size 2 21 , 
Ax = 1, At — 0.1, D — a = 0.25, initially, = everywhere). 

ized by the well-known 0, critical exponents of the 
GV class. A direct integration of Eq.(|5J fulfills all these 
expectations: a space-time plot of reveals the typi- 
cal branching-annihilating dynamics of kinks (Fig. |3t) 
schematized in Fig. At criticality, the density of in- 
terfaces decays with an exponent 9 — 0.28(1) fully com- 
patible with the expected value 0.286 (Fig. Eb). 

In three dimensions, preliminary results indicate that 
the scenario observed in d = 2 remains valid: for large 
b > 0, separate Ising and DP transitions are found, while 
for (essentially) negative b values, a unique, direct transi- 
tion occurs. This last transition, however, does not pos- 
sess the marginal character of its two-dimensional coun- 
terpart: indeed it is "fully first-order", in the sense that 
the density of interfaces now also jumps discontinuously 
to zero. We note that this is in agreement with d = 2 
being the upper critical dimension for the GV class. 

We now summarize and comment on our results. The 
Langevin equation (0) does account for all the known 
phenomenology about order-disorder phase transitions in 
the presence of two symmetric absorbing states. In d = 1, 
such a description puts an end to the discussion, in micro- 
scopic models, about the relevance of "parity conserva- 
tion" . We note en passant that in the absorbing phase of 
the GV transition our equation becomes an effective con- 
tinuous description of the annihilation process 2 A — > 0, 
as testified by the l/Vt decay in Fig. Oh- In d = 2, we ob- 
tained clear evidence that a full GV line exists, but splits 
into an Ising and a DP line at some parameter value b* 
(inset of Fig. i n qualitative agreement with mean- 
field predictions. In d = 3, and more generally above 
the upper critical dimension of the GV class, the same 
scenario is found but for the fact that the GV transition 
is now fully first-order. 



It is not clear to us whether Eq. J5J can be derived 
from first-principles, especially given the importance, in 
microscopic models exhibiting the GV transition, of the 
annihilation reaction 2 A — > 0, for which this is reput- 
edly impossible [2l| . Nevertheless, both the generating 
functional associated to Eq. and the one rigorously 
obtainable from the master equation for the correspond- 
ing reaction-diffusion processes enjoy the same symmetry 
and feature similar characteristic invariants. Thus ana- 
lytical studies of our proposal appear promising, be it 
either within renormalisation-group perturbative calcu- 
lations (maybe akin to those of 221), or via the nonper- 
turbative approach put forward in [5j . First results under 
the latter auspices are very encouraging 0. 
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